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ABSTRACT 

We present the study of the emission properties of the quasar PG1211+143, which belongs 
to the class of Narrow Line Seyfert 1 galaxies. On the basis of observational data analyzed by 
us and collected from the literature, we study the temporal and spectral variability of the source 
in the optical/UV/X-ray bands and we propose a model that explains the spectrum emitted in 
this broad energy range. In this model, the intrinsic emission originating in the warm skin of the 
accretion disk is responsible for the spectral component that is dominant in the softest X-ray 
range. The shape of reflected spectrum as well as Fe K line detected in hard X-rays require the 
reflecting medium to be mildly ionized (£ ~ 500) . We identify this reflector with the warm skin of 
the disk and we show that the heating of the skin is consistent with the classical aP to t prescription, 
while aPgas option is at least two orders of magnitude too low to provide the required heating. 
We find that the mass of the central black hole is relatively small [Mbh ~ 10 7 — 10 8 Af Q ), which 
is consistent with the Broad Line Region mapping results and characteristic for NLS1 class. 

Subject headings: accretion, accretion disks - quasars: individual (PG 1211+143) - X-rays:galaxies - 
X-rays:spectra - X-rays:variability 



1. Introduction 

The quasar PG1211+143 (z = 0.0809; Marziani 
et al. 1996) is one of the prototypes of AGN with 
particularly strong hard optical spectrum and pro- 
found steep soft X-ray emission (Elvis, Wilkes, & 
Tananbaum 1985). It has been frequently argued 
that the optical/UV and soft X-ray emission form 
a single Big Blue Bump component which extends 
across the unobserved XUV band (Bechtold et al. 
1987), and strongly dominates the bolometric lu- 
minosity. 

However, the nature of the soft X-ray emission 
- in particular, the spectral component that ex- 
hibits a very soft spectrum and often is seen in 
many quasars and Seyfert 1 galaxies only below 



~ 1 keV - is still not determined unambiguously. 
There are arguments both in favor and against 
the origin of this emission component being the 
tail of the optical/UV Big Blue Bump. This is 
because any interpretation relies on extrapolating 
the spectrum over a decade in frequency through 
the unobserved XUV band. 

The strong arguments in favor of the Big Blue 
Bump (BBB) extending up to the soft X-ray range 
are statistical in their nature. There is a strong, 
almost linear correlation between the luminosity 
of an object at 2500 A and the soft X-ray lumi- 
nosity in the ROSAT spectral band both in a large 
sample of Seyfert 1 galaxies (Walter & Fink 1993) 
and in radio quiet quasars (Yuan et al. 1998) sug- 
gesting a relatively uniform spectral shape of the 
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BBB. Spectral analysis of both high rcdshift and 
low redshift quasars (Laor et al. 1997) allowed for 
a "shrinking" of the unobserved band and such 
spectra appear to cover smoothly the entire opti- 
cal/UV/soft X-ray band. 

Arguments against the soft X-ray emission be- 
ing a direct tail of the optical / UV component are 
also mostly statistical. It is rather surprising from 
the theoretical point of view why the extension of 
the BBB should be the same in all objects having 
a range of black hole masses and luminosities. In- 
stead, it was suggestive of some atomic processes 
(e.g. Czerny & Zycki 1994). Also the normalized 
variability amplitude in the UV band is usually 
smaller than in the X-ray band (see e.g. Ulrich et 
al. 1997). A close look at the UV spectra of some 
sources frequently suggests a slight decline on vF v 
diagram in the far UV (Malkan 1983; Magdziarz 
et al. 1998 for NGC 5548). This docs not really 
allow for this component to extend far enough into 
XUV range to make it observable in the soft X- 
ray band, even if a certain amount of unsaturated 
Comptonization is taken into account. 

Detailed analysis of the broad band spectra 
of Seyfert 1 galaxy NGC 5548 (Magdziarz et al. 
1998) led to a conclusion that the optical/UV 
spectrum is well modeled by an accretion disk, 
and the hard X-rays are reproduced by standard 
thermal Comptonization. However, an additional 
component is needed to model the soft X-ray 
source: the most viable candidate here is another, 
optically thick Comptonizing medium possibly as- 
sociated with the transition between the disk and 
the hot plasma. This component would contain 
a major part of the source bolometric luminosity. 
Such a model also supports the scenario includ- 
ing a strong and extended Big Blue Bump with 
Comptonization playing a major role in formation 
of this component. 

In this paper we study in detail the quasar 
PG121 1+143. Our goal is to investigate (i) if a 
similar model also applies to this source, which 
is representative both for quasars and for Narrow 
Line Seyfert 1 galaxies, (ii) if this model is unique 
and (iii) what is the physical basis and possible 
observational consequences of such a phenomeno- 
logical model. 

The choice of object was motivated by the avail- 
ability of considerable amount of observational 
data available, which span the last two decades. 



The object can be considered as a typical, al- 
though rather faint quasar, with \F\ at 5100 
A about 5 x 10 44 erg s" 1 (Kaspi et al. 2000) 
for the Hubble constant H Q = 75 km s _1 Mpc -1 . 
The broad band a ox index measured between 2500 
A and 2 keV is 1.16 (Corbin & Boroson 1996) 
while for a typical quasar, it is 1.5 (Yuan et al. 

1998) . The soft X-ray slope measured in the 
ROSAT band is - 3.1 (Walter & Fink 1993) while 
typical quasar value is 2.58 (Yuan et al. 1998). Its 
optical emission lines are narrow, with the FWHM 
of Hp line is 1832 ± 81 km s" 1 (Kaspi et al. 2000; 
see also Corbin & Boroson 1996; Wilkes et al. 

1999) implying a classification of this object as 
a Narrow Line Seyfert 1 galaxy. 

In Section 2 we describe the observational data 
that we analyzed, and how we supplement these 
with the data available from the literature. In Sec- 
tion 3 we give the results of the variability and 
spectral studies of PG121 1+143. We discuss the 
physical processes that may lead to the observed 
behavior of the source in Section 4 and summarize 
the results in Section 5. 

2. Observations and data analysis 

First spectral X-ray data for PG1211+143 were 
obtained from the EINSTEIN observatory (see Ta- 
ble 1). Strong soft X-ray emission seen in EIN- 
STEIN was described by Bechtold et al. (1987), 
and it was subsequently reanalyzed by Elvis et al. 
(1991). We determine the 0.2 - 2 keV and 2 - 
10 keV intrinsic fluxes (i.e. corrected for Galactic 
reddening) in the second observation from the fits 
provided by Bechtold et al. (1987). The quality 
of the first observation was much lower, so we de- 
termined the soft and hard fluxes in this data by 
rescaling the Bechtold et al. (1987) results with 
the IPC count rate ratios 0.695 between the two 
observations, as given by Elvis et al. (1991). 

Subsequent observations were performed by 
EXOSAT on three occasions (Elvis et al. 1991; 
Saxton et al. 1993). We determine the 2 - 10 keV 
intrinsic fluxes through the spectral fitting of pub- 
licly available ME data to a simple broken power 
law model with fixed galactic absorption. The 0.2 
- 2 keV flux was calculated from the results of 
spectral fitting to LE data presented in Saxton et 
al. (1993). 

On one occasion the source was observed by 
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Table 1 
x-ray observations 



Instr. 


Date 


F(0.2 - 2) a 


F(2 - 10) a 


EINSTEIN 


1979Dec05 


52.1 


7.3 


EINSTEIN 


1980Decll 


75.0 


10.5 


EXOSAT 


1985Junl3 


59.7 


9.9 


EXOSAT 


1986Jan28 


41.2 


6.4 


EXOSAT 


1986Jan09 


18.9 


5.8 


GINGA 


1985May06 




8.5 


ROSATa 


1991Decl7 


16.1 




ROSATb 


1992Junl7 


21.9 




ROSATc 


1992Jun28 


30.6 




ROSATd 


1992Jun28 


34.2 




ROSATe 


1993Jun02 


20.5 




ROSATf 


1993Jun25 


19.0 




ASCA 


1993Jun03 


8.5 


3.7 


RXTE 


1997Augl6 




12.9 



a Fluxes without Galactic absorption in units of 
10~ 12 erg s^ 1 cmT 2 



GINGA. We take the 2 - 10 keV flux from Lawson 
& Turner (1997). 

PG121 1+143 has been observed six times by 
the ROSAT satellite with the PSPC detector. One 
of the ROSAT pointings was simultaneous with 
ASCA observations. The ROSAT and ASCA 
spectral data were reduced using the standard 
software package, where for ASCA, we used data 
from all four detectors (SIS0/1 and GIS2/3) fit- 
ted simultaneously. Preliminary results based on 
these data were presented by Yaqoob ct al. (1994) 
but we repeat their analysis, using the updated, 
current version of the ASCA response matrices. 
The ASCA observation was sufficiently long (88 
ksec) to attempt a study of the light curve within 
this single data set, so we could also determine the 
variability on short time scales. 

The Rossi X-ray Timing Explorer (RXTE) ob- 
servations of PG1211+143 are new, conducted as 
a part of optical / X-ray monitoring program. 
The object was monitored roughly every 5 days 
for 6 months; the results of this program will 
be reported elsewhere (Netzer, Madejski, Kaspi, 
et al., in preparation). Here, we use the spec- 



trum derived from the summed RXTE Propor- 
tional Counter Array (PCA) data analyzed in a 
standard manner (including the background sub- 
traction), using the ftool script rex. The source 
was detected to at least 12 keV, beyond which 
the signal - to - noise ratio was becoming very 
low. There was no significant detection with the 
RXTE's High Energy X-ray Timing Experiment 
(HEXTE). 

All analysis of the spectral X-ray data was done 
using the XSPEC software package, version 10.0 
(Arnaud 1996). We adopt a redshift 0.0809 to 
this source after Marziani et al. (1996) instead of 
0.085 used in most previous papers (e.g. Bechtold 
et al. 1987; Yaqoob et al. 1994) since the new 
determination seems to be convincing. Whenever 
we fix the Galactic column density, we adopt the 
value of Nh = 2.83 x 10 20 cnr 2 after Elvis, Wilkes, 
& Lockman (1989) although another value (Nh = 
2.1 x 10 20 cm~ 2 ) was used by Lawson & Turner 
(1997); we note that the exact value of Nh makes 
essentially no difference in our analysis. 

The optical/UV spectrum of this source was de- 
termined by Bechtold et al. (1987). We use the 
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broad band data of Elvis et al. (1994) and we com- 
pare them with the HST observations of Bechtold 
et al. (2000), performed on April 13 and 16, 1991. 

3. Results 

3.1. Time variability 

3.1.1. Short timescales in the X-ray band 

We study the short timescale variability in the 
hard and soft X-ray bands on the basis of the 
longest sequence from ASCA (June 3 - 4, 1993, 
8.7 x 10 4 s), and accompanying ROSAT (June 2 
- 4, 1993, 4 x 10 4 s). We analyzed the ASCA 
lightcurve binned into 128 s bins for two SIS and 
two GIS instruments separately, in the entire en- 
ergy band. As an example, we show the GIS-2 
lightcurve in Figure 1. 

The source is clearly variable in both energy 
bands. The rms amplitude in ASCA and ROSAT 
bands are given in Table 2. Such values are compa- 
rable to those measured in X-rays in other AGN, 
e.g. 16% in MCG-6-30-15 on timescales shorter 
than 10 6 s (Nowak & Chiang 2000). The count 
rate in ASCA is dominated by 1 - 2 keV energy 
band, where both the contribution of the soft X- 
ray excess and hard X-ray power law are impor- 
tant. We therefore compared the rms value de- 
rived above with the rms value measured from 
the light curve for photons above 2 keV. We used 
only the data from GIS instruments as those data 
are more reliable for time-series analysis than SIS 
data for PG 1211+143. This is because the PG 
1211+143 ASCA observation was conducted in a 
4-CCD mode of the SIS instrument, and the source 
image fell close to the gap between the CCD chips. 
Since the satellite attitude can change somewhat 
during the observation, some unknown fraction of 
photons may fall in the space between the chips, 
resulting in an apparent but erroneous variabil- 
ity. Both GIS2 light curves (whole energy band, 
and above 2 keV) were rebinned into 512 s bins 
to increase statistics; the two observations yield 
rms values of 0.095 and 0.171, respectively. This 
suggests that the observed variations at the time 
scales of minutes-hours are predominantly due to 
the hard X-ray power law. This is confirmed by 
the fact that rms value measured in the ROSAT 
band is much smaller during the 40 ksec observa- 
tion than in the ASCA data. 



The data from ASCA and ROSAT satellites 
are exactly simultaneous; nonetheless, the pres- 
ence of gaps due to the Earth occultation of the 
source makes an attempt to use the Fourier anal- 
ysis rather difficult, although there are methods 
which allow to overcome this problem for high 
quality observations (e.g. Done et al. 1992a). 
However, even the available ASCA data are not 
of sufficient quality to perform this kind of anal- 
ysis. Therefore, here we adopt a complementary 
method based on a nonlinear prediction, which is 
described in Czerny & Lchto (1997). We analyzed 
the light curve binned into 128 s bins for two SIS 
and two GIS instruments separately, in the entire 
energy band. Such analysis reveals that the shape 
of the power density spectrum (hereafter PDS) is 
basically described as a power law, with the slope 
1.5. Errors in determination of those values are 
substantial, and systematic errors dominate, since 
the existence of the signal in the data is marginal 
(the correlation coefficient itself is 0.62 and 0.67, 
i.e. below 0.8, which would be necessary for reli- 
able results). 

Such a value of the PDS slope is typical for 
Seyfert galaxies at high frequencies. Yaqoob 
(1997) obtained the high frequency slope of MCG- 
6-30-15 of 1.4, Edelson & Nandra (1999) deter- 
mined the slope to be 1.74 in NGC 3516 at fre- 
quencies above ~ 10 -6 Hz, and Nowak & Chi- 
ang (2000) as well as Chiang et al. (2000) give 
the slope of 1.0 at intermediate frequencies and a 
slope of 2.0 at high frequencies in MCG-6-30-15 
and NGC 5548, respectively. 

We conclude that the variability properties of 
PG121 1+143 are rather typical of a Seyfert galaxy. 
Unfortunately, neither quasars nor Narrow Line 
Seyfert 1 galaxies were studied with equal accu- 
racy as normal Seyfert 1 galaxies, so it is not 
known whether there are any systematic differ- 
ences in variability apart from scaling with mass 
(e.g. Leighly 1999; Piasecki et al. 2000). Results 
for PG1211+143 may indicate that the variability 
pattern is universal. 

3.1.2. Trends of long-term variability in the X- 
ray band 

Available archival data for PG1211+143 in- 
clude six ROSAT pointings as well as historical 
observations from EINSTEIN, EXOSAT, GINGA, 
and ASCA; with those, we can also study the 
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Table 2 

Variability properties of PG1211+143 



Timescale 


F(0.2 - 2) Q 


rms 


F(2 - 10) a 


rms 


3.3 x 10 4 


17.9 









8.7 x 10 4 






3.96 


0.09 


1.7 x 10 5 


18.2 


0.06 






4.8 x 10 7 


23.7 


0.25 






4.3 x 10 8 


36.5 


0.62 


7.44 


0.30 


5.6 x 10 8 






8.12 


0.34 



a Fluxes in units of 10 12 erg s 1 cm 2 



0.2 




2xl0 4 



4xl0 4 6x10" 8x10" 
Time (s) 



Fig. 1.— The ASCA GIS-2 lightcurve in 2 - 10 
keV energy band. 



-0.2 



-0.4 




0.2 0.4 0.6 0.8 1 
log (AT/128s) 



Fig. 2. — The dependence of the decay of correla- 
tions with the time step in the ASCA lightcurve. 
First three points were used to determine the slope 
(solid line) which directly translates into the slope 
of the power density spectrum. 
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long time variability patterns of the source. The 
ROSAT observations cover a period of 1.5 years, 
allowing for a reliable measurement with the same 
instrument. We determine the observed fluxes in 
0.2-2 keV band assuming a simple spectral model 
of a broken power law, with the second (hard X- 
ray) power law photon index fixed at 2.0 and the 
Galactic absorption fixed at the assumed value as 
given above. The ROSAT sequences c and d ob- 
tained on the same day were both fitted separately 
and combined (see Table 3). Since the errors in 
separate fits are large, we use only the combined 
set in further analysis. All flux measurements are 
summarized in Table 1. 

The variability during the period of ROSAT 
observations is moderate, roughly by a factor 2. 
More precisely, the rms value measured on the 
timescale of 4.8 x 10 7 s is only 0.25, which is not 
much higher than the rms variability on short time 
scales, within a single ASCA data set. 

However, at the timescales of years the trend 
reverses, and the soft X-ray variability amplitude 
becomes higher than the hard X-ray variability on 
the same time scale. Analyzing the entire period 
of twenty years of observations, we measure the 
rms in 0.2 - 2 keV band which is higher than in 2 
- 10 keV band (see Table 2). It is consistent with 
the trends shown by the Big Blue Bump in other 
AGN (e.g. Clavel et al. 1992 for NGC 5548 and 
Clavel, Wamsteker & Glass 1989 for Fairall 9). 

The most dramatic change happened during the 
end of 1985/bcginning of 1986 but this period is 
covered by only 3 EXOSAT observations (see Fig- 
ure 3). Large changes in the soft X-ray emission 
are accompanied by the changes in the hard X— 
ray flux but not with comparable amplitude (see 
Figure 4). 

Since the analysis of ASCA data allows us to 
estimate the slope of the PDS at high frequencies, 
we can reconstruct the total PDS from the values 
of rms on time scales of 8 x 10 4 s and 5 x 10 s s. We 
assume the slope of 1.5 at high frequencies and a 
slope 0.0 below certain frequency v a . Using the in- 
tegral conditions connecting the PDS and rms at a 
given frequency (see e.g. Czerny, Schwarzenberg- 
Czerny, & Loska 1999) we obtain the flattening 
frequency v Q = 8.4 x 10~ 7 Hz, corresponding to a 
timescale of 1.2 x 10 6 s. This is again comparable 
to the results for most Seyfert galaxies: the flat 
part of the PDS starts at timescales ~ 10 6 s for 



MCG-6-30-15 (Nowak & Chiang 2000), - 10 7 s 
for NGC 5548 (Chiang et al. 2000) and - 10 8 s 
for NGC 3516 (Edelson & Nandra 1999). 

3.1.3. Black hole mass determination from vari- 
ability 

X-ray variability time scales have been success- 
fully used to estimate the mass of the central black 
hole in active galaxies (e.g. Barr & Mushotzky 
1986; Wandel & Mushotzky 1986), although a 
more robust method is the measurement of fea- 
tures in the PDS (see, e.g., Edelson & Nandra 
1999). Convenient and powerful method based on 
the normalization of the high frequency tail of the 
power density spectrum was recently developed by 
Hayashida et al. (1998) and Hayashida (2000). 
Here wc use a modified version of this method (Pi- 
asecki et al. 2000) which is based on comparison of 
the shape of the PDS with independent mass de- 
termination for the Seyfert 1 NGC 5548, an AGN; 
for the purpose of the mass estimate, this is likely 
to be more appropriate than scaling the normaliza- 
tion of Cyg X-l, an accreting binary system, used 
in Hayashida et al. 1998. The data discussed in 
previous section are not of sufficient quality to per- 
form Fourier analysis. However, having the rms 
determined for a few time scales, we can recon- 
struct the power density spectrum assuming its 
shape in advance and using the integral connec- 
tion between the power density spectrum and rms 
(see e.g. Czerny et al. 1999). 

Here we assume that the PDS is fiat below cer- 
tain frequency and has a slope 1.5 above it, as de- 
termined in Section 3.1.1. As data points for this 
method, we use the GIS rms measurement for the 
8.7 x 10 4 s time scale and multi-instrument mea- 
surements on the 4.3 x 10 8 s time scales. Extension 
of the measurements afforded by the RXTE obser- 
vations only confirms that the break frequency is 
at timescales below ~ 10 8 s; if this timescale is 
above ~ 10 8 s, the rms would saturate. With this, 
we determine the break frequency to be at 2 x 10~ 7 
Hz, corresponding to timescales of 5 x 10 6 s. The 
error on this value is large, since the uncertainty 
of the GIS rms measurement is already a factor 5. 

This simple approach gives the value of the 
black hole mass logM = 7.0 ± 0.7 [M©], if the er- 
ror in determination of the rms in GIS data (0.3) 
is taken into account. The mass of PG1211+143 
is thus definitely rather small for a quasar, and 
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Fig. 3. — The time dependence of the soft 0.2 - 2 
keV flux from instruments: EINSTEIN, EXOSAT, 
ASCA/ROSAT and ROSAT. 
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Fig. 4. — The soft 0.2-2 keV flux as a function 
of 2 - 10 KeV hard X-ray flux from various in- 
struments: EINSTEIN (filled squares), EXOSAT 
(filled circles), ASCA/ROSAT (open circle), and 
ROSAT (open squares). 



more typical for Seyfert galaxies. High luminos- 
ity of the object has to be related to a relatively 
high accretion rate, close to the Eddington value. 
Small mass is consistent with the measurements 
of Kaspi et al. (2000) based on the properties of 



the Broad Line Region: M 
L/L E dd = 0.29. 
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3.1.4- Spectral variability in the soft X-ray band 

While the overall X-ray spectral variability 
as discussed above can be interpreted MS ("I rela- 
tive change of the intensity of the soft vs. the 
hard component, we also aim to determine if the 
soft component alone of the total spectrum of 
PG121 1+143 has a variable shape. To this end, 
we fit the ROSAT data with a spectral model 
consisting of a cold absorber and a broken power 
law. The results are shown in Table 3. We also 
searched for the presence of an ionized absorber 
in this source by adding an ABSORI XSPEC model 
in the case of the data set /, which shows the 
worst x 2 . The fit was only marginally improved 
(Ax 2 = 7.0, with 34 d.o.f. instead of 36) and the 
best fit value of the column of the ionized gas was 
5.5 x 10 21 cm~ 2 (between 1.5 and 13.0), with a 
temperature kT ~ 1 x 10 6 K and £ ~ 20, where 
both T and £ arc only weakly constrained. 

The data sets clearly show the variability of the 
flux of the soft component, as discussed in the pre- 
vious section. Those variations, however, are not 
accompanied by significant changes in the spectral 
shape. There is a slight trend in the variations of 
the model value of Nh with the total flux. How- 
ever, no such trend is visible in the photon index. 
Formal analysis of the Njj — F v relation shows 
that the two parameters are best represented by a 
linear fit shown in Fig. 5 but a constant line rep- 
resentation is also acceptable (Nh — 3.44 ± 0.17, 
X 2 = 0.60 for 6 d.o.f.) since the errors in determi- 
nation of Nh are very large (see Table 3). Similar 
analysis of the r — F v trend confirms the visual 
impression that the correlation between the two 
quantities is insignificant - and a constant line fit 
gives T = 3.21 ± 0.07 (x 2 = 2.6 for 6 d.o.f.). With 
this, we conclude that the shape of the soft X-ray 
component does not vary, within the accuracy of 
the available data. 
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Table 3 

Fits of broken powerlaw model to ROSAT data 



Data N H l\ E break T 2 X 2 /d.o.f. Flux(0.2 - 2) a 



a 


°-^°-0.35 


3.0918:2 


1.8o_q 50 


2.0 


39.0/36 


16.13 


b 


4fi +0.46 
°.*°-0.42 


q qq + 0.18 

J -°°-0.17 


> 2.0 6 


2.0 


39.6/36 


21.94 


c 


q qc+0.65 
o - oo -0.60 


q 99+0.24 
°- zz -0.23 


>2.0 


2.0 


30.3/36 


30.61 


d 


q qq + 0.79 

o - oo -0.70 


q 1 7 +0.30 
J - i ' -0.28 


1 41 +°° 


2.0 


48.2/36 


34.25 


cd 


q 9 n +0.44 
d - zu -0.44 


3.n±8:iS 


i . D '-0.31 


2.0 


48.5/36 


30.92 


e 


GC+0.43 
°- oo -0.40 


3.3818:3 


>2.0 


2.0 


53.2/36 


20.53 


f 


cq+0.49 
°- ou -0.33 


3.1818:11 


1-781S°35 


2.0 


62.3/36 


18.99 



a Fluxes without Galactic absorption 
b Energy break beyond the ROSAT band 



Table 4 

Fits to combined ROSAT data with variable normalization 



Model Pari" Par2 b Par3 c Par4 d Par5 e X 2 /d.o.f. 



A bknpower 


c 4 +0.18 


q 9 r+0.07 

°- z '- ) -0.06 


> 1.75 


2.0 




306/192 


B power 


ro+0.18 
°- OJ -0.17 


q 9C+0.07 
°- z '- ) -0.06 








306/191 


C bbody+power 


1 01 +0.13 
i - ol -0.12 


n 107 +0.003 
U - 1U( -0.003 




2.0 




331/191 


D comptt 


2.83 


< 0.020 


3Q+ - 07 

u ' OB -0 .06 




10.4+?-8 


294/191 


E comptt+power 


2.83 


< 0.017 


24+ ' 11 


2.0 


14.61*;$ 


290/190 



a Galactic absorption in units 10 20 cm 2 

b Soft spectral index (models A and B) or the soft photon temperature in keV (models C 
and D) 

c Brcak energy (model A) or the electron temperature (model D) 
d Hard X spectral index (models A, C and E) 
e Optical depth (models D and E) 
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3.1.5. Optical/UV variability 



3.6 - 



3.4 



3.2 




10 20 30 40 

F (0.2 - 2 keV)/10" 12 erg/s/cm z 



Fig. 5. — The relation between the fitted value of 
N H and the soft 0.2 - 2 kcV flux in the ROSAT 
data. Dotted line represents the best fit constant 
value. Dashed line represents the best linear fit. 
Errors in Nh are very large (~ 0.4; see Table 3). 
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Fig. 6. — The relation between the value of photon 
index V and the soft 0.2 - 2 kcV flux in the ROSAT 
data. Dotted line represents the best fit constant 
value. Dashed line represents the best linear fit. 
Errors in T are large (~ 0.2; see Table 3). 



PG121 1+143 was extensively monitored in the 
optical band in the period 1966 - 1982, but no 
variability was found within the accuracy of ob- 
servations (Barbieri & Romano 1984). Bcchtold et 
al. (1987) found the UV flux from IUE measure- 
ments in 1981 and 1982 that was systematically 
lower than the trend suggested by the optical data, 
although with a similar slope. Wit this, they grey- 
shifted the UV data by +0.27 in logarithmic scale 
to match the the optical data. More recent HST 
observations (Dobrzycki, private communication) 
give exactly the same UV flux as the grey-shifted 
value of Bechtold et al. (1987) (-25.68 of Dobrzy- 
cki at 1387 A vs. -25.69 of Bechtold at 1370-1424 
A in the same logarithmic units), making it dif- 
ficult to distinguish the true variability from an 
instrumental effect. The question was finally set- 
tled by the monitoring performed over a 7 1/2 year 
time span by Kaspi et al. (2000). New optical data 
clearly showed the variability by a factor of two, 
and the normalized variability amplitude at 5100 
A was determined to be 16.2%. 

3.2. Spectral fitting of the broad-band 
data 

Here, we perform spectral fits using simple phc- 
nomenological models as well as more realistic, 
physical models to individual data sets, and then 
make an attempt to fit the broad-band data si- 
multaneously. We note that despite the changes in 
the normalization of each component, the shape of 
each spectral component varies very weakly with 
time. With this, we attempt to describe this shape 
more accurately by using the combined data for 
each component. We assume all model parameters 
are equal in all data sets apart from normalization. 

3.2.1. ROSAT data 

Firstly, we analyze the combined ROSAT PSPC 
data, which we fit simultaneously using the same 
spectral model with variable normalization. We 
fit the following phenomcnological models: broken 
power law with a variable break energy but with 
fixed hard photon index of 2.0 (model A); a single 
power law (model B); and black body with an un- 
derlying power law of photon index 2.0 (model C). 
We also fit more physical models, namely: thermal 
Comptonization model (comptt) of Titarchuk 
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(1994) parameterized by the soft photon tempera- 
ture, hot plasma temperature and its optical depth 
(model D). The best fitting results are presented 
in Table 4. 

The last model clearly provides the best de- 
scription of all the ROSAT data, although the 
X 2 is still quite high. The soft photon temper- 
ature has only an upper limit of 0.034 keV, or 
4 x 10 5 K. Such temperature is consistent with 
the effective temperature of a standard accretion 
disk around a black hole with the mass inferred 
by us above. The hot plasma temperature is rela- 
tively low, kT e ~ 0.4 keV or 5 x 10 6 K, much lower 
than the plasma responsible for the formation of 
the hard X-ray emission, and the optical depth of 
the plasma is very high, r ~ 10. The additional 
power law tail added in the model does not im- 
prove this fit substantially, since the ROSAT data 
cannot constrain the spectral slope in hard X-rays. 

We also tried to improve the quality of fits 
to the spectrum by including the in the model 
the OVIII emission line. The model parameters 
(model E in table 4) remained basically unchanged 
and the line energy was fitted to 653 eV. However, 
the fit was not improved (x 2 /d.o.f. =287/189), and 
the line equivalent width could not be constrained 
(EW = 32± 2 ° eV). 

3.2.2. Combined ROSAT/ASCA data - model 
without reflection 

We apply our best model inferred above to 
ASCA/ROSAT data, which include all ROSAT 
pointings; in the fits we use data from all four 
ASCA detectors, as described previously. We sup- 
plement the COMPTT component with an under- 
lying hard X-ray power law component character- 
ized by its photon index, and assume fixed Galac- 
tic absorption with N H = 2.8 x 10 20 cm~ 2 . The 
results are given in Table 5, and the spectra are 
displayed in Figure 7. 

The fit to the data in the entire ROSAT/ASCA 
energy band is acceptable, and the Comptoniz- 
ing plasma parameters are only slightly modified 
such that the Comptonizning plasma temperature 
is now lower, and the optical depth increased by 
a factor of 2. However, the basic picture re- 
mains qualitatively unchanged. The hard X-ray 
index is close to typical values for Seyfert galaxies 
and somewhat smaller than typically measured in 



quasars. In this fit the reflection component was 
neglected. 

The results for the combined ROSAT/ASCA 
data differ slightly from the results obtained on the 
basis of the ASCA data alone by Vaughan et al. 
(1999). The model used by Vaughan et al. is much 
simpler (power law + black body). Nevertheless, 
it is clear that the slope of the power law in our fit 
is flatter than in theirs (1.88 versus 2.07) and the 
soft X-ray excess is relatively stronger - the ratio 
of the luminosity of the soft X-ray excess to that 
of the hard power law in the 0.6 - 10 keV ASCA 
band is 0.28 in our fits and 0.17 in theirs. This 
difference is probably caused by our inclusion of 
the ROSAT data and can be attributed to residual 
calibration uncertainties between ROSAT PSPC 
and ASCA, resulting in different spectral slopes 
(see Iwasawa et al. 1999, describing the analysis 
of nearly-simultaneous ROSAT PSPC and ASCA 
data for NGC 5548). 

3.2.3. ASCA data - model with reflection 

The spectral fitting of the ASCA data alone 
were presented by Vaughan et al. (1999) for the 
case of a simple model consisting of a power law, a 
black body and a Gaussian modeling the iron Ka 
line. While the description of the soft X-ray excess 
with a simple black body component is justified 
in the ASCA data - these data are not sensitive 
to the shape of the soft excess - the presence of 
X-ray reflection in the data was not adequately 
studied by these authors, although the detection 
of the iron K a line in the ASCA data (Yaqoob et 
al. 1994; Vaughan et al. 1999) clearly shows its 
importance. 

Current models for the structure of AGN inter- 
pret the presence of the iron K a line resulting from 
the irradiation of the accretion disk surface by 
hard X-ray emission. Such a fluorescent line is ac- 
companied by the presence of a reflection compo- 
nent which is not an independent spectral feature 
(see, e.g., George & Fabian 1991). Therefore, we 
considered the properties of the entire reflection 
component, i.e. reflection hump plus the iron line 
(Zycki, Done, & Smith 1997; 1998). The repro- 
cessed spectrum is parameterized by the ioniza- 
tion parameter £ = Lx/n e r 2 (Done et al. 1992b) 
and the reflection amplitude R = Cl/2w, where fl 
is the solid angle subtended by the reprocessor as 
viewed from the X-ray source. We applied two 
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Table 5 

Results of simultaneous spectral fitting to 
ASCA and ROSAT data 



r 


K1 soft 


kT a 




1.88 ±0.05 


< 0.02 


0.15 ±0.01 


21 ±2 1.18 



a Temperature in keV 



Note. — The model consists of the Comptonized 
black body, power law and Galactic absorption (fixed 
value) . 



versions of this model: one that does not contain 
the effect of relativistic smearing and another that 
includes this effect. The second model formally 
contains additional free parameters, but some of 
these parameters of the model were always held 
fixed. In our case, the power law index of the ra- 
dial distribution of the incident flux was fixed at 
-3.0, the entire disk surface was assumed to be a 
reflector (i.e. Ri n = 6R g , R out = 10 3 i? g , where 
Rg = GM/c 2 ), and the inclination angle was as- 
sumed to be 30°. The effect of these studies is 
shown in Table 6. 

Formal minimum \ 2 is f° r the relativistically 
smeared neutral reflection, but the normalization 
of the reflection component is surprisingly high. 
The assumption about the model adopted for the 
soft X-ray excess was not essential - fitting the 
COMPTT model instead of a black body led to the 
same values of relativistic reflection. If the amount 
of reflection is fixed at 1.0 the reflection is possi- 
bly ionized (£ = 4001^°) but £ is not well con- 
strained. If the reflection is assumed to be neutral, 
the amplitude of reflection is high (R = 2.8^ ' 9 ). 
The quality of the fits is the same as before, and 
the \ 2 f° r the second option is lower by 3.0. 

3.2.4- RXTE data - reflection component 

In order to put better constraints on the re- 
flection parameters, we used the RXTE data, 
which nominally extend up to w 20 keV. While 
PG1211+143 in its flux state, as observed with the 
RXTE, is a relatively faint source to be studied in 
great detail by the RXTE's PCA, a non-imaging 
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channel energy (keV) 



7. — The model of the combined 
fly ASCA data (see Table 5). 
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Table 6 

Results of fitting the reflection component to ASCA data. 



Refl T kT b % £ R X 2 /d.o.f. 



non relat 


2.03 ± 0.05 


101 ± 0.12 




442.4/448 


non rclat 


2.15 ± 0.07 


93 ± 0.14 




431.1/449 


non relat 


2.11 ± 0.07 


95 ± 0.14 500 


0.45 +0 3 6 


436.4/4 49 


relat 


2.18 ± 0.08 


91 ± 0.15 6 601 66™ 


1.0 


429.2/ 449 


relat 


2.28 ± 0.07 


85 ± 0.14 


3.611:? 


425.1/449 



a Temperature in eV 

Note. — The model consists of black body, power law, Galactic absorption 
(fixed value) and reflection component with an iron line. 



instrument, the data still put some constraints on 
the crucial 8-20 keV band where the Compton 
reflection is important. We adopted the compos- 
ite model consisting of a primary soft component 
(described as a Comptonized black body), power 
law continuum, and the Compton - reflected spec- 
trum with the iron Ka line included (Zycki, Done 
& Smith 1997; 1998). The reprocessed spectrum 
is parameterized by £, R — Q/2ir, as above. The 
outer disk radius, as before, was fixed at 10 4 R g 
and the power law index describing the incident 
flux as a function of radius was fixed again at the 
value -3.0. However, the inclination angle and the 
inner radius of the reflecting disk surface were now 
allowed to vary. We only studied the relativisti- 
cally smeared model since it is the most physical 
one. The iron line energy, again, is not a free pa- 
rameter, but it depends on the ionization state 
of the reflecting medium. Its profile is connected 
with the assumed geometrical parameters: inner 
radius of accretion disc and viewing angle. 

The resulting best fit parameters are given in 
Table 7 and presented in Figure 8. New results did 
not differ significantly from the ROSAT/ASCA re- 
sults regarding the soft part of the spectrum, al- 
though the optical depth of the Comptonizning 
medium dropped somewhat, and the hard X-ray 
power law became softer. However, they provide 
an interesting insight into the geometry. The am- 
plitude of reflection is consistent with 1.0, and the 
reflecting area extends down to the marginally sta- 



> 
> 




2 5 
channel energy (keV) 



Fig. 8. — The model fit (continuous histogram) to 
the combined non-simultaneous ASCA/XTE data 
(see Table 7). Dashed histogram shows the con- 
tribution of the soft COMPTT component, dotted- 
dashed histogram - the reflection component and 
dots mark the hard X-ray power law. 
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ble orbit which means that the disk is probably not 
disrupted. However, its surface (at least in the 
inner parts dominating the reflection) is consider- 
ably ionized (£ <~ 500), and the neutral reflection 
is excluded at a 90% confidence level. For such a 
case, the resulting iron line energy corresponds to 
the domination of Fe XXII ions and its rest frame 
energy is E ~ 6.7 keV, while the line equivalent 
width is EW ^150 cV. This supports the results 
of Yaqoob et al. (1994), who found a marginal evi- 
dence for an ionized iron Ka line fitting the ASCA 
data. The inclination angle of the disk is low, i.e. 
we see the disk essentially face on. 

3.2.5. Optical/ UV/ X-ray data 

Finally, we include the optical/UV data into 
our modeling. For that purpose, we add a new 
model component, which essentially has a form 
of standard accretion disk of Shakura & Sunyaev 
(1973) around a non-rotating black hole. We as- 
sume that the disk, whether it is or it is not cov- 
ered by the ionized skin, radiates locally as a black 
body. It is certainly a simplifying assumption (see 
e.g. Merloni, Fabian & Ross 2000) but more accu- 
rate models critically depend on the assumed de- 
tails of the disk vertical structure (e.g. Nayakshin, 
Kazanas, & Kallman 2000; Madej & Rozaiiska 
2000b). Such assumptions are beyond the scope 
of the present paper. 

Fits of the reflected component suggest that 
the disk extends down to the marginally stable or- 
bit. However, the reflection component indicates 
a highly ionized hot disk surface, and thus we as- 
sume that the emission of the standard cold disk 
is only detectable from radii larger than certain 
Ri on . The gravitational energy dissipated above 
Ri on is directly re-emitted by the disk, while the 
energy dissipated below Ri on has to provide the 
energy source for hard X-ray emission and Comp- 
tonizing medium. 

The disk parameters are therefore: the black 
hole mass M, the mass accretion rate M, and 
Rion- We assume a face-on view of the accretion 
disk, and the distance to the source calculated as- 
suming H a = 75 km s _1 Mpc _1 . 

In the upper panel of the Figure 9 we plot 
the broad band spectrum of PG1211+143, con- 
structed for the case of the maximum value of tem- 
perature of the seed photons undergoing Comp- 



tonization. The solid line marks the sum of the 
cold accretion disc spectrum and the model fit- 
ted to the X-ray data that consists of the Comp- 
tonized black body (dashed line), power law func- 
tion and Galactic absorption (which is neglected 
in the plot). 

In this case there is a clear gap between the 
UV and soft X-ray band, which would be seen as 
a local minimum in the continuum around 2 x 10~ 2 
keV. The bolomctric luminosity measured between 
10 000 A and 100 keV (assuming an extrapola- 
tion of the measured hard X-ray power law) is 
1.5 x 10 45 erg s _1 . The mass of the black hole is 
2.2 x 10 8 M Q , the accretion rate is 0.4M© yr _1 , 
and R ion — 28R g . L/LEdd ratio in this solution 
is 0.05. About 50 % of the energy is dissipated in 
the black body component of the disk. 

In the lower panel of Figure 9 we show the 
broad band spectrum determined under the as- 
sumption that cold disk itself provides the seed 
photons for Comptonization, i.e. the soft pho- 
ton temperature is given by the effective temper- 
ature of the inner disk. Such a spectrum forms 
a continuous Big Blue Bump extending from the 
optical/UV to soft X-ray band. The bolomet- 
ric luminosity of the source is in this case higher, 
1.9 x 10 45 erg s _1 . The mass of the black hole is 
1.5 x 10 8 M Q , the accretion rate is 0.54 M Q yr _1 , 
and Ri on — A2R g . L/LEdd ratio in this solution is 
equal to 0.10. About 30 % of the energy is directly 
emitted by a black body outer disk, while 60 % is 
in the Comptonized component and ~ 10 % in the 
hard X-ray power law. 

Both fits presented above reproduce the broad 
band spectral data well. We note that all inter- 
mediate solutions are also possible and at a level 
of pure data analysis, we cannot reject any of the 
above possibilities. However, a composite spec- 
trum of quasars given by Laor et al. (1997) and 
the broad band spectrum of RE J1034+396 (e.g. 
Puchnarewicz et al. 2000) show clearly the exis- 
tence of a single broad Big Blue Bump extending 
from optical/UV to soft X-ray band. Therefore 
using the analogy with these data we would favor 
our second solution. 

The mass of the central black hole determined 
from the spectral analysis is somewhat higher than 
the value determined from the variability analysis 
(see Sect. 3.1.3). However, all quoted values are 
rather small for a quasar and comparable to the 
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Table 7 

Results of simultaneous fitting to ASCA and RXTE data 



r 


Kl soft 


kT a 


T 


R 


c 


cos i 






2.18 ±0.03 


< 0.03 


0.13 ±0.05 


12 ±5 


85+ - 65 


r nf) +600 

ouu_ 450 


0.9 ±0.1 


6.02^-0 


1.04 



a Tcmpcrature in keV 
b Inner radius in R g 



Note. — The model consists of the Comptonized black body, power law, Galactic absorption (fixed value) 
and reflection component with an iron line. 



masses of the central black holes in Seyfert 1 galax- 
ies such as NGC 5548, despite significantly higher 
bolomctric luminosity. 

4. Discussion 

4.1. Determination of the black hole mass 

The mass and luminosity estimates presented 
in this paper and those available in the literature 
render PG1211+143asa relatively low-luminosity, 
low-mass object. Quasars constituting the sample 
of Laor et al. (1997) and Zhcn ct al. (1997) have 
black hole masses on average about 1.4 x 10 9 M Q 
but they have also higher luminosities, about 10 46 
erg s _1 . An accurate measurement, however, is 
quite difficult to obtain. Our method based on 
the X-ray variability in the form of PDS measure- 
ment (see Sect. 3.1.3) gives log(M) = 7.0 ± 0.7, 
while the spectral fit indicates log(M) = 8.17. We 
did not give the formal error of this second mea- 
surement since the systematic error certainly dom- 
inates any statistical error involved. Nonetheless, 
the error in the value of the mass determined from 
the spectral fitting can be roughly estimated under 
the assumption that the outer part of the flow is 
well described by a standard Keplerian disk. The 
fit is based on two relationships. The first one is 
the asymptotic shape of the disk emission in opti- 
cal band 

F v oc ^^{MMfl^co^Hlf-^, (1) 

where cos i is the inclination angle of the disk, H Q 
is the adopted value of the Hubble constant and / 



is the spectral correction given by the color tem- 
perature to the effective temperature ratio. The 
second constraint comes from the bolometric lu- 
minosity of the system 

F bol cxMr)Hl (2) 

where r\ is the efficiency of accretion. In our fits 
we adopted H a = 75 km s -1 Mpc -1 , cosi = 1.0, 
/ = 1 and 77 = 1/16. However, if those values 
in reality are different from the assumed ones, the 
derived mass will scale as 

M = 1.5 x 10 8 / 2 r7(cosi)- 1 (i/ o /75)" 1 . (3) 

This means that the departure from a black 
body approximation and slight inclination of the 
system will tend to increase the mass, although 
not considerably. The color to effective tempera- 
ture ratio was determined to be ~ 1.7 (Shimura 
& Takahara 1995; Sobczak et al. 1999). Merloni, 
Fabian, & Ross (2000) give even larger values for 
Galactic sources, but in the case of AGN, efficient 
bound-free transitions in heavy elements reduce / 
significantly (Rozaiiska 2000, private communica- 
tion). Therefore the upper limit for / seems to 
be about 1.7, and the most probable inclination 
is 30°, giving together an increase in the derived 
mass by a factor of 3. On the other hand, if there 
is a strong outflow close to the black hole, it will 
effectively decrease the standard efficiency of ac- 
cretion, thus decreasing the mass. Strong outflow 
may give a factor of two, if it is comparable ener- 
getically to the outflows in radio-loud objects. On 
the grounds of this analysis, we can estimate an 
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kT .,,,.=0.018 keV 
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Fig. 9. — The broad band spectrum of 
PG121 1+143 in the IR - X range. In the lower 
panel we show the model calculated assuming 
that the accretion disk provides seed photons for 
Comptonization by the hot plasma. In the up- 
per panel we show the model which assumes that 
the soft photons are emitted in the hotter plasma. 
The solid triangles are data points from Elvis et 
al. (1994). The UV spectrum with emission lines 
are the HST data (Bechtold et al. 2000) and 
the crosses are the X-ray data from ASCA and 
ROSAT. The Comptonized disk model was fitted 
to the X-ray data. In both cases, absorption in 
soft X-rays is neglected in the plot, and for the 
modeling purposes, the X-ray data were corrected 
for absorption. Optical/UV spectrum was mod- 
eled with a standard disk. Thick solid line is the 
sum of all model components. 



error in our determination from the spectral anal- 
ysis as log(M) = 8.17t°: J. The lower limit is still 
not consistent within a formal error with the de- 
termination based on variability, but in that case 
it is difficult to estimate the possible systematic 
error of the method. As for the determination of 
Kaspi et al. (2000), the formal errors are small: 
4-05 ti'^i x 10 7 M Q , but recent analysis of the sys- 
tematic errors of methods based on BLR proper- 
ties indicates a possible factor of three error in 
both directions, due to technical problems with 
sampling as well as a possible (and quite proba- 
ble) flattening of the BLR (Krolik 2000). 

All this implies that taking the systematic er- 
rors into account, all measurements are roughly 
consistent with each other, even though they are 
highly inaccurate. We notice, however, that if we 
allowed for a maximally rotating black hole as a 
central object, our mass based on the spectral fits 
would be higher by a factor of 7. This result, how- 
ever, would be rather difficult to reconcile with 
other measurements. 

4.2. Luminosity state 

PG121 1+143 is both a quasar, and a Nar- 
row Line Seyfert 1 galaxy. Temporal and spec- 
tral analysis presented in Section 3 show that the 
mass of the central black hole is rather low (~ 
5x 10 7 — 2x 10 8 M Q ) for a source of such optical/X- 
ray bolometric luminosity (~ 1.5 — 1.9 x 10 45 erg 
s ). In a sample of Kaspi et al. (2000) it is on 
the border between normal Seyfert galaxies, NLS1 
and quasars, having mass typical for a Seyfertl but 
narrow Hp line and larger luminosity than nor- 
mal Seyferts. An interesting comparison can be 
made with Fairall 9 galaxy having broad Hp: when 
PG1211+143 is faint and Fairall 9 is exception- 
ally bright (Clavel, Wamsteker & Glass 1989) their 
luminosities are nearly equal. For those cases, 
the luminosity in Eddington units L/LEdd of PG 
1211+143 remains only marginally larger than for 
Fairall 9 (we adopted the black hole mass of Firall 
9 to be twice that of PG121 1+143, after Kaspi et 
al. 2000). 

The exact value of L/LEdd is not determined 
accurately since it depends on the black hole mass 
and bolometric luminosity determination. Our 
L/LEdd ratio determined from spectral fits (~ 
0.05—0.1) is somewhat lower than the value of 0.29 
given by Kaspi et al. (2000) due to larger black 
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hole mass and lower bolomctric luminosity than 
resulted from BLR study and the use of monochro- 
matic luminosity at 5100 A. The source is strongly 
variable and our fits (as well as results of Kaspi et 
al. 2000) are mostly representative for a relatively 
faint state. 

This result supports the view expressed by 
many authors (e.g. Pounds et al. 1995; Pounds 
& Vaughan 2000) that Narrow Line Seyfert 1 
galaxies are objects accreting at a relatively high 
L/L<Edd ratio. The analysis of the Compton re- 
flection implies that the inclination to the accre- 
tion disk surrounding the nucleus is low, which is 
contrary to the scenario where NLS1 galaxies are 
viewed at a high inclination and their variability is 
caused by obscuration events (Brandt & Gallagher 
2000; see also Boiler 2000). 

Weak hard X-ray producing corona (as com- 
pared to the large luminosity of the Blue Bump) 
and large L/ LEdd ratio are consistent with the pre- 
diction of the two-temperature corona model (Ja- 
niuk & Czerny 2000), and the relatively low pho- 
ton index even in this source (r = 2.18) suggests 
that models with strong advection cannot repro- 
duce the observed spectrum well (Janiuk, Zycki, & 
Czerny 2000). However, a non-accreting magnetic 
corona model, perhaps including some dumpiness, 
may well be viable, as argued by Poutanen & 
Fabian (1999). 

4.3. Reflecting and Comptonizing medium 

Our spectral analysis requires the existence of 
a number of various plasma zones: (1) very hot 
plasma, responsible for the hard X-ray emission, 

(2) moderately hot, moderately dense plasma re- 
sponsible for Comptonization, (3) moderately hot 
plasma responsible for reflection and (4) a cold, 
dense disk. It is natural to pose the question 
whether the components (2) and (3) are actually 
the same medium but parameterized differently in 
the spectral fits. 

Medium (2) has the temperature ~ 0.15 keV, or 
10 6 K, and a large optical depth, r ~ 10. Medium 

(3) has the ionization parameter £ ~ 500. Tem- 
peratures intermediate between the black body 
temperature and the Inverse Compton tempera- 
ture arc difficult to achieve by radiatively heated 
plasma, since there is a rapid transition between 
the two solutions (Krolik, McKee & Tarter 1981). 



This transition happens at S w 10 (see, e.g., Zycki 
& Czerny 1994), and the relationship between the 
two ionization parameters is £ = 5 x 10~ 5 TS. 
Plasma in zones (2) and (3) should be in such a 
transition state. Assuming S = 10 and £ = 500 
from plasma zone (3) we obtain the temperature 
consistent with the temperature of plasma zone 
(2). This would imply that the same medium 
which Comptonizes the disk photons coming from 
below also reflects the hard X-ray photons illumi- 
nating it from above. 

Such interpretation clearly simplifies the geom- 
etry of the flow. The accreting stream consists 
simply of the outer black body disk and an inner 
disk which develops very optically thick (r <~ 10), 
warm skin (T ~ 10 6 K). A small fraction of en- 
ergy is released in the form of hard X-ray emission 
above this skin. A possible schematic view of the 
flow is shown in Figure 10. 

Similar values of the optical depth and the tem- 
perature of the "warm" plasma were derived for 
NGC 5548 by Magdziarz et al. (1998). However, 
in this source the magnitude of reflection was lower 
and it was not possible to locate geometrically the 
Comptonizing medium. Magdziarz et al. (1998) 
actually suggested that the radial transition zone 
from a classical disk to an inner ADAF-likc flow is 
the most probable location. Another example of 
the presence of the warm skin was inferred from 
the XMM observations of PKS 0558-504 (O'Brien 
et al. 2000). Interestingly, numerous examples of 
a Comptonizing medium that modifies the shape 
of the spectrum in soft X-ray band were inferred 
from data for Galactic sources: for those, even be- 
fore the effects of Comptonization, this soft com- 
ponent extends to higher energies as compared 
with AGN due to the higher disk temperature, and 
therefore can be more easily observed. Examples 
include Nova Muscae 1991 (Zycki et al. 1998), Cyg 
X-l (Gierlihski et al. 1999), GRS 1915+105 (Vilhu 
& Nevalainen 1998), GS 2000+25, XTE J1550- 
564, GRO J1655-40 (Zycki et al., in preparation), 
and the problem has been discussed in a review of 
Zycki (2000). 

Our results seem to indicate that, at least 
in sources like PG1211+143, this warm medium 
forms a horizontal layer. Its optical depth is only 
t ~ 10, which is rather low in comparison to 
the total optical depth of the classical disk. Still, 
about 40 - 60% of the total gravitational energy 
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dissipated by the flow is released in this zone, while 
the remaining <~ 50 — 30% is released in the body 
of the disk. Roughly 10% of energy is emitted in 
hard X-rays and thus presumably released in the 
hot corona. However, this quantity is difficult to 
estimate due to the unobserved high energy cutoff. 

4.4. Formation of the warm skin 

Application of the thermal instability of the ir- 
radiated hot plasma by Krolik et al. (1981) to the 
vertical structure of the disk clearly shows that 
a warm skin, or a corona at the Inverse Comp- 
ton temperature, forms above the cold disk. How- 
ever, the details of the sharp transition between 
the disk and the Inverse Compton - heated corona 
are not clear (e.g. Begelman, McKee, & Shields 
1983; Rozahska & Czerny 1996; 2000; Nayakshin, 
Kazanas, & Kallman 2000; Nayakshin 2000). 

Current models imply that the optical depth of 
this transition zone of intermediate temperature is 
generally very low and the total optical depth of 
the heated skin is never large. Results of Nayak- 
shin ct al. (2000) and Zycki & Rozahska (2000) 
allow at most for r ~ 1.5 even if the X-ray radi- 
ation is released in compact regions leading to an 
enhancement of the irradiation under a hot spot 
by as much as two orders of magnitude. This is 
simply related to the fact thatthe X-ray heating 
does not penetrate deeply enough to support the 
high temperature of the transition zone. 

Since the Comptonized component contains a 
significant fraction of the bolomctric luminosity 
(~ 40 — 60 %), much larger than the fraction con- 
tained in the hard X-ray part (up to ~ 10 %) 
the Comptonizing medium cannot be heated by 
hard X-rays. For this reason, the discrepancy be- 
tween the results quoted above and the large op- 
tical depth determined from our fits is not sur- 
prising. The dissipation must occur within this 
zone and we can attempt to compare the required 
amount of heating with most popular parametric 
prescription for the viscous dissipation of gravita- 
tional energy. 

The properties of the reflection and Comp- 
tonization by the warm skin allow us to estimate 
the number density and the temperature of this 
medium at the optical depth of r <~ 1.0. The tem- 
perature T warm ~ 10 6 K results directly from the 
COMPTT fit to the data. The density can be deter- 



mined from the ionization parameter of reflection 
£ = 500, assuming that the typical radius under 
discussion is 20R g and the bolometric luminosity 
in hard X-ray band is 10 44 erg s _1 : 

n= ^^l 12 cm -3 (4) 
r ? 

or equivalently, p ~ 10~ 12 g cm~ 3 . This medium 
is clearly Compton-cooled, since the efficiency of 
bremsstrahlung for such parameters, namely 

k B = 6.6 x l0 20 p 2 T 1/2 ~ 0.006 erg s _1 cm _3 (5) 

is two orders of magnitude lower than efficiency of 
Comptonization, which is 

AkT 

A c = F tot K es p ? <~ 0.2 erg s _1 cm~ 3 , (6) 

where the value of F tot ~ 10 15 ergs _1 cm -2 was es- 
timated from the total bolomctric luminosity 10 45 
erg s _1 and the distance of 20R g . 

This large amount of cooling has to be compen- 
sated by the heating. The most popular param- 
eterization of the heating in accretion flow is the 
a viscosity prescription introduced by Shakura & 
Sunyaev (1973). Two variants are most popular: 
aPtot and aPgas- In the first case the amount 
of heating (in the context of a Kcplcrian flow) is 
given by 

Qptot = l-SaP-totClK ~ 0.3a erg s _1 cm~ 3 , (7) 

where we approximated Ptot by P ra d ~ Ftot / c; the 
Kcplcrian angular velocity Qk was estimated at 
20R g for a 10 8 M Q black hole mass. In the second 
case 

Qp a as = ^aP gas fl K ~ 3a x 10~ 3 erg s _1 cm _3 , 

(8) 

where the density and the temperature were taken 
from the previous estimates. 

We see that the aP gas variant of heating is by 
two orders of magnitude too small to compensate 
for the cooling. However, the aP tot variant is of 
the right order of magnitude, provided the viscos- 
ity coefficient a is large. It does not explain the 
physical mechanism behind the heating but clearly 
indicates that aPtot scaling is a good phenomeno- 
logical description of the process. 
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4.5. Consequences of the aP tot heating 

There are two major consequences of the scaling 
favored by our analysis for the dissipation: disk 
instability and surface temperature gradient. 

The disk models are known to be unstable un- 
der such scaling if the radiation pressure domi- 
nates, as it is the case in disks in active galaxies 
(Pringle, Rees, & Pacholczyk 1973; Ligthman & 
Eardley 1974). The disk is supposed to display 
strong outbursts in timescales of years, with about 
two orders of magnitude difference between the 
low and high luminosity states, as it is suggested 
by direct scaling of computations of the time evo- 
lution of Galactic sources (e.g. Honma ct al. 1991; 
Szuszkiewicz & Miller 1998; Janiuk ct al. 2000). 

The variations of the Big Blue Bump on these 
timescales are observed in active galaxies in gen- 
eral (see Sect. 3.1.2), and the maximum variability 
(peak on PDS diagram) has not been determined 
yet for the far UV and soft X-ray band. This may 
be because such a peak is possibly at even longer 
timescales than can be inferred from observations 
covering a span of only 20 years. Therefore, we 
cannot yet determine the viscosity parameter a 
governing the timescales of outbursts. One aspect 
of the current data which may be in conflict with 
the above scenario is that the observed variations 
of the source luminosity do not appear as a frac- 
tion of a global outburst with the predicted large 
amplitude. There arc, however, various effects 
which may lead to the difference between the mod- 
eled outbursts and the observed time behavior. 
The models studied so far did not include a warm 
skin as that discussed in Sect. 4.4. Moreover, any 
strong outflow is known to suppress the outburst 
amplitude, if the outflow is coupled to the tempo- 
rary local accretion rate (Janiuk et al. 2000), and 
any irregularities in this outflow lead immediately 
to much more complex lightcurve than simple pre- 
dictions would indicate (Nayakshin, Rappaport, 
& Melia 2000). In addition, if the viscosity pa- 
rameter depends on the amount of the radiation 
pressure, the amplitude of the outburst can be re- 
duced. Testing such scenarios is more promising 
in case of Galactic sources because of the shorter 
timescalesrequired . 

Another aspect of the scenario that may present 
a conflict with observations is connected with the 
broadly discussed issue of the lack of Lyman edge 



in the observed spectra of AGN while the pres- 
ence of such an edge is predicted by most models 
of accretion disk atmospheres (for a review, see 
Koratkar & Blaes 1999). Those studies, however, 
never included such a strong heating close to the 
surface as well as the accompanying Comptoniza- 
tion. Strong heating leads to a more isothermal 
atmosphere, as already argued by Laor & Netzer 
(1989), and Comptonization is also quite efficient 
in removing the edge, as stressed by Czerny & 
Zbyszewska (1991). Some of the recent papers al- 
ready made a step towards reducing the predicted 
edge (Rozariska et al. 1999; Madej & Rozanska 
2000a; Hubeny & Hubeny 1998) but most likely 
the issue cannot be fully resolved until heating is 
properly included in computations of the disk at- 
mosphere. 

4.6. Future observational tests 

High resolution spectroscopy in soft X-ray band 
available now through Chandra and XMM data 
may provide an independent test of the physical 
conditions in the reflecting/Comptonizing warm 
skin. This region is not completely ionized so it 
is a source of intense emission lines, mostly of the 
hydrogen-like Lya lines from CNO as well as Si, 
Mg, S and Ne. The presence of such lines has been 
possibly detected in the MCG-6-30-15 spectrum 
taken with the XMM (Branduardi-Raymont et al. 
2000). 

Computations of the irradiated slab predict the 
existence of such lines (Ross, Fabian, & Young 
1999). Dumont et al. (2000) performed the radia- 
tive transfer in an irradiated slab under constant 
pressure, with the spectral slope of the incident 
radiation r = 2.0 and the high energy cut-off of 
100 keV. In this case the equivalent width of the 
Lya CVI, NVII and OVIII lines were equal to 1.6 
cV, 2.3 eV and 7.4 eV for the ionization parameter 
£ (at the surface) of 300, and 1.3 cV, 2.8 eV and 
11.5 cV for £ = 1000. Those values are quite simi- 
lar to the EW measured in the X-ray spectrum of 
MCG-6-30-15 (Branduardi-Raymont ct al. 2000). 
However, since in those considerations the slab was 
only radiatively heated from the outside, the high 
temperature zone was relatively thin (Thompson 
depth of 0.1 and 0.3, correspondingly). The op- 
tical depth of the Comptonizing medium deter- 
mined by the data was on the order of 10. With 
this, it is difficult to predict if the additional hcat- 
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ing in the interior of the warm skin interior may 
change these results. 

On the other hand, the theoretical computa- 
tions performed for a constant density irradiated 
medium, with harder incident flux (T = 1.9, ex- 
tending up to 280 keV) indicated much larger EW, 
of order of 50 eV for £ = 1000 (Czerny & Dumont 
1998) which suggests that EW of the soft X-ray 
lines is sensitive to the structure of the warm skin. 
Clearly, more detailed theoretical studies together 
with observational constraints on the line intensity 
and its profiles may help to resolve the problem 
of the formation of such warm skin. The mod- 
els should include both the effect of Comptoniza- 
tion and kinematical/gravitational broadening of 
the line profiles, as both processes lead to devel- 
opment of red wings in emission lines, and may be 
difficult to differentiate. 

5. Conclusions 

The Narrow Line quasar PG1211+143 is con- 
sidered to be a very good candidate to study the 
broad band emission and variability properties due 
to a rich sample of observational data available for 
this object. In this work, we considered the phys- 
ical mechanisms which may be responsible for the 
optical, UV and X-ray emission and variability of 
this source, and we found that the following model 
can give a satisfactory explanation of the observa- 
tions: 

• The central object powering the source is a 
supermassive black hole with a mass M = 
10 7 - 2 x 10 8 M Q and the luminosity to Ed- 
dington luminosity ratio is in the range of 5 
- 40 %. Despite large errors both quantities 
are consistent with the explanation of the 
phenomenon of Narrow Line Seyfert galaxies 
based on a relatively low mass of the black 
hole and a high accretion rate. 

• The emitted radiation originates in the cold 
accretion disk responsible for the optical/UV 
spectrum. Inner disk region develops a 
warm, optically thick skin, which produces 
the profound "soft excess" observed in the 
X-ray spectrum and is well modeled with 
a Comptonized black body emission. How- 
ever, from the spectral analysis alone, it 
cannot be determined whether the source of 



seed photons for Comptonization is the cold 
accretion disk or a somewhat hotter plasma, 
and both explanations are possible. 

• The warm disk skin is also responsible for the 
reflection of hard X-rays, which are emitted 
in the accretion disk corona. The latter has 
the form of active flares rather than forming 
a continuous medium. 

• Both accretion disk and its warm skin are 
heated by the viscous dissipation process, 
which can be parameterized by otPtot pre- 
scription. This may lead to the source lu- 
minosity variations resulting from the pres- 
ence of regions where the radiation pressure 
is dominant. 

• Measurements of the soft X-ray emission 
lines possible in the Chandra and XMM data 
and their comparison with theoretical mod- 
els may help to test the warm skin structure. 
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Fig. 10. — The geometry of the accretion flow in 
PG1211+143 consistent with the spectral model 
presented in Fig. 9. The optical flux is emitted by 
the cold accretion disk (T ~ 10 4 K). The disk is 
the source of seed photons for the hot Comptoniz- 
ing cloud (T ~ 10 6 K, r ~ 20), which extends 
below the transition radius Ri on . The hard X-ray 
flux is emitted by the hot flare region (T ~ 10 9 K) 
and are partially reflected by the cloud (£ ~ 500, 
fl/2ir ~ 1). 
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